An analysis of over-relaxation in a kinetic approximation of systems of conservation laws

The over-relaxation approach is an alternative to the Jin–Xin relaxation method in order to apply the equilibrium source term in a more precise way. This is also a key ingredient of the lattice Boltzmann method for achieving second-order accuracy. In this work, we provide an analysis of the over-relaxation kinetic scheme. We compute its equivalent equation, which is particularly useful for devising stable boundary conditions for the hidden kinetic variables.     

Liouville theorem for weighted p-Lichnerowicz equation on smooth metric measure space

In this paper, let $(M^n,g,d\mu )$ be n-dimensional noncompact metric measure space which satisfies Poincaré inequality with some Ricci curvature condition. We obtain a Liouville theorem for positive weak solutions to weighted p-Lichnerowicz equation

${△}_{p,f}v+c{v}^{\sigma }=0,$

where $c\le 0,m>n\ge 1,1<p<\frac{m?1+\sqrt{(m?1)(m+3)}}{2},\sigma \le p?1$ are real constants.     

A probabilistic approach to spectral analysis of growth-fragmentation equations

The growth-fragmentation equation describes a system of growing and dividing particles, and arises in models of cell division, protein polymerisation and even telecommunications protocols. Several important questions about the equation concern the asymptotic behaviour of solutions at large times: at what rate do they converge to zero or infinity, and what does the asymptotic profile of the solutions look like? Does the rescaled solution converge to its asymptotic profile at an exponential speed? These questions have traditionally been studied using analytic techniques such as entropy methods or splitting of operators. In this work, we present a probabilistic approach: we use a Feynman–Kac formula to relate the solution of the growth-fragmentation equation to the semigroup of a Markov process, and characterise the rate of decay or growth in terms of this process. We then identify the Malthus exponent and the asymptotic profile in terms of a related Markov process, and give a spectral interpretation in terms of the growth-fragmentation operator and its dual.     

Mass distribution and skewness for passive scalar transport in pipes with polygonal and smooth cross sections

We extend our previous results characterizing the loading properties of a diffusing passive scalar advected by a laminar shear flow in ducts and channels to more general cross‐sectional shapes, including regular polygons and smoothed corner ducts originating from deformations of ellipses. For the case of the triangle and localized, cross‐wise uniform initial distributions, short‐time skewness is calculated exactly to be positive, while long‐time asymptotics shows it to be negative. Monte Carlo simulations confirm these predictions, and document the timescale for sign change. The equilateral triangle appears to be the only regular polygon with this property—all others possess positive skewness at all times. Alternatively, closed‐form flow solutions can be constructed for smooth deformations of ellipses, and illustrate how both nonzero short‐time skewness and the possibility of multiple sign switching in time is unrelated to domain corners. Exact conditions relating the median and the skewness to the mean are developed which guarantee when the sign for the skewness implies front (more mass to the right of the mean) or back (more mass to the left of the mean) “loading” properties of the evolving tracer distribution along the pipe. Short‐ and long‐time asymptotics confirm this condition, and Monte Carlo simulations verify this at all times. The simulations are also used to examine the role of corners and boundaries on the distribution for short‐time evolution of point source , as opposed to cross‐wise uniform, initial data.     

Local linear convergence analysis of Primal–Dual splitting methods

In this paper, we study the local linear convergence properties of a versatile class of Primal–Dual splitting methods for minimizing composite non-smooth convex optimization problems. Under the assumption that the non-smooth components of the problem are partly smooth relative to smooth manifolds, we present a unified local convergence analysis framework for these methods. More precisely, in our framework, we first show that (i) the sequences generated by Primal–Dual splitting methods identify a pair of primal and dual smooth manifolds in a finite number of iterations, and then (ii) enter a local linear convergence regime, which is characterized based on the structure of the underlying active smooth manifolds. We also show how our results for Primal–Dual splitting can be specialized to cover existing ones on Forward–Backward splitting and Douglas–Rachford splitting/ADMM (alternating direction methods of multipliers). Moreover, based on these obtained local convergence analysis result, several practical acceleration techniques are discussed. To exemplify the usefulness of the obtained result, we consider several concrete numerical experiments arising from fields including signal/image processing, inverse problems and machine learning. The demonstration not only verifies the local linear convergence behaviour of Primal–Dual splitting methods, but also the insights on how to accelerate them in practice.     

A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.     

Evolution of weak noise and regular waves on dissipative shock fronts described by the Burgers model

The interaction of weak noise and regular signals with a shock wave having a finite width is studied in the framework of the Burgers equation model. The temporal realization of the random process located behind the front approaches it at supersonic speed. In the process of moving to the front, the intensity of noise decreases and the correlation time increases. In the central region of the shock front, noise reveals non-trivial behaviour. For large acoustic Reynolds numbers the average intensity can increase and reach a maximum value at a definite distance. The behaviour of statistical characteristics is studied using linearized Burgers equation with variable coefficients reducible to an autonomous equation. This model allows one to take into account not only the finite width of the front, but the attenuation and diverse character of initial profiles and spectra as well. Analytical solutions of this equation are derived. Interaction of regular signals of complex shape with the front is studied by numerical methods. Some illustrative examples of ongoing processes are given. Among possible applications, the controlling the spectra of signals, in particular, noise suppression by irradiating it with shocks or sawtooth waves can be mentioned.     

三维位势问题的梯度边界积分方程的新解法

三维位势问题的边界元分析中,关于坐标变量的边界位势梯度的计算是一个困难的问题. 已有一些方法着手解决这个问题,然而,这些方法需要复杂的理论推导和大量的数值计算. 本文提出求解一般边界位势梯度边界积分方程的辅助边值问题法. 该方法构造了与原边界值问题具有相同解域的辅助边值问题,该辅助边值问题具有已知解,因此通过求解此辅助边值问题,可获得梯度边界积分方程对应的系统矩阵,然后将此系统矩阵应用于求解原边值问题,求解过程非常简单,只需求解一个线性系统即可获得原边值问题的解. 值得注意的是,在求解原边值问题时,不再需要重新计算系统矩阵,因此辅助边值问题法的效率并不很差. 辅助边值问题法避免了强奇异积分的计算,具有数学理论简单、程序设计容易、计算精度高等优点,为坐标变量梯度边界积分方程的求解提供了一个新的途径. 3个标准的数值算例验证了方法的有效性.      

A survey of surface tension,molar volume and density for Sn–Ag–Cu–Bi–Sb quinary alloys as lead-free solders

ABSTRACT

By applying the geometric models and the theoretical equation, the surface tension, the molar volume and the density were studied. The empirical calculations were carried out in temperature range 623?K?≤?T?≤?1123?K. Only few thermophysical properties were estimated for eight quinary alloys: Sn3.55Ag0.5Cu3Bi3Sb, Sn3.48Ag0.5Cu3Bi5Sb, Sn3.48 Ag0.5Cu5Bi3Sb, Sn3.40 Ag0.5Cu5Bi5Sb, Sn3.53Ag1Cu3Bi3Sb, Sn3.46Ag1Cu3Bi5Sb, Sn3.46Ag1Cu5Bi3Sb, Sn3.38Ag1Cu5Bi5Sb. The results show that surface tension and density have a linear appearance for all temperatures. We have also studied the influence of the composition and temperature in the studied alloys. The obtained theoretical results are compared with the experimental ones and with the conventional Pb–Sn welds.